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ABSTRACT 

Context. HgMn Chemically Peculiar stars are among the quietest stars of the main-sequence. However, according to theoretical 
predictions, these stars could have pulsations related to the very strong overabundances of iron peak elements, which are produced by 
atomic diffusion in upper layers. Such pulsations have never been detected from ground based observations. 

Aims. Our aim is to search for signatures of pulsations in HgMn stars using the high quality lightcurves provided by the CoRoT 
satellite. 

Methods. We identified three faint stars (V > 12), from VLT-GIRAFFE multiobject spectrograph survey in a field which was planned 
for observation by CoRoT. They present the typical characteristics of HgMn stars. They were observed by the CoRoT satellite during 
the long run (131 days) which started from the 24th of October 2007, with the exoplanets CCD's (Additional Programme). In the 
present work, we present the analysis of the ground based spectra of these three stars and the analysis of the corresponding CoRoT 
lightcurves. 

Results. Two of these three HgMn candidates show low amplitude (less than 1.6 mmag) periodic variations (4.3 and 2.53 days 
respectively, with harmonics) which are compatible with periods predicted by theoretical models. 

Key words. Stars: chemically peculiar - Stars: oscillations (including pulsations) - Techniques: photometric 



1. Introduction 

HgMn stars are known to form a group among the Chemically 
Peculiar stars of the main-sequence (Preston, 1974). Their 
photosphere is characterized by strong abundance anoma- 
lies of many metals (see for instance Castelli&Hubrig 2004). 
Overabundances with respect to solar ones can be about two 
orders of magnitude and even larger. For instance, Mn can be 
enhanced by more a factor of 200, and Hg by a factor up to lO*". 
These stars are found in the range of effective temperature be- 
tween about lOOOOK and 16000K (Pop I, B7-A0 IV- V), they are 
slow rotators, often binaries, and no magnetic fields has been 
yet detected (Shorlin et al. 2002, Wade et al. 2006). These stars, 
which are supposed to have no convective motions in their atmo- 
sphere, which do not show any indication of activity, are possibly 
among the quietest stars of the main sequence. 

The trends of abundance peculiarities observed in HgMn 
stars are usually well understood in the framework of the atomic 
diffusion (first proposed by Michaud, 1970). As stated by the 
classical model, diffusion process builds up abundance stratifica- 
tions in the stable atmosphere (and envelope) of HgMn stars, ac- 
cording to the balance between radiative accelerations and grav- 

Send offprint requests to: G. Alecian 

* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 



ity. Theoretical predictions by Alecian & Michaud, (1981) for 
Mn were in excellent agreement with observations (Smith & 
Dworetsky, 1993). However, modelling is not yet able to deal 
with the large variety of peculiarities observed in these stars. On 
another hand, detailed studies of stratification build up remain to 
be done, because of the heaviness of numerical computations for 
optically thin media (see Alecian & Stiff, 2007). Modelling of 
atomic diffusion in stellar interior (optically thick case) is more 
advanced (Richard et al. 2001). Even if these studies of stellar in- 
teriors cannot address the problem of atmospheric peculiarities, 
they allow to precise the evolution of these stars. A recent work 
by Turcotte & Richard (2003) predicts that, due to iron accumu- 
lation in the upper layers of their envelopes, these stars should at 
least be slow pulsators (as the Slow Pulsating B stars). Presently, 
there is no observational evidence of pulsation in HgMn stars, 
and CoRoT satellite offers an exceptional opportunity to con- 
strain the models. 

A detailed description of the CoRoT mission can be found in 
the CoRoTbook (Baglin, 2006). Beside the two CCD's devoted 
to asteroseismologjQ the exoplanet CCDs photometrically mon- 
itor thousands of stars in the magnitude range 10.5 < niR < 16 . 
This monitoring is designed to detect planet transits but it has 
also a high interest for studies of stellar photometric variations. 



' Each of these CCD's, can measure photometric variation of 5 bright 
stars (oti/ < 9.5) 
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Stellar science outside the CoRoT core programme have been 
then allowed to access the data of exoplanet channel through the 
Additional Programme procedure (Weiss, 2006). 

In Sect. |2] we present our strategy to find, through ground 
based observations (discussed in Sect. [3]), faint HgMn stars in 
fields observed by CoRoT. In Sect. ID we present the analysis 
of the CoRoT data for the targets we identified as HgMn can- 
didates. Finally, we describe theoretical models for HgMn stars, 
followed by a discussion in view of our observational results. 



2. Observational strategy 

According to our knowledge of the solar neighborhood, one 
can estimate that about 8% of main-sequence AV-BV stars are 
HgMn stars. Therefore, applying this rate to the density of main- 
sequence AV-BV stars predicted by the "Besangon" model of 
stellar population synthesis of the Galaxy (Robin et al. 2003), 
in the anticenter direction, only about 0.3% of Pop I stars of 
any type may be HgMn. Then, one can expect to have about 35 
HgMn stars in a field containing 1 1400 stars, which is approxi- 
mately the number of stars monitored by exoplanet CCD's. The 
problem is that there is not any available catalog of such faint 
HgMn stars, since the MK classification with the dispersion suf- 
ficient to recognize the HgMn signature was caiTied out only for 
relatively bright stars. 

We have then decided to realize prior to the launch of 
CoRolQ a survey in fields pre-selected for exoplanets search 
in the direction of the galactic anticenter, and using VLT- 
FLAMES/GIRAFFE (ESQ) multiobject spectrograph. 

The VLT-GIRAFFE 25' field towards the anticenter direc- 
tion contains about 490 IV-V stars 11.5 < my < 16, according 
to the "Besangon" model. The instrument can provide spectra of 
1 10-120 objects in each field, with good enough resolution and 
S/N to identify abundance peculiarities (see next section). If stars 
were selected at random, the probability to find HgMn stars were 
extremely weak. To overcome this difficulty, we used a prelimi- 
nary version of the Exo-Dat catalog (Deleuil et al. 2008), where 
the stellar types in our fields were roughly determined in view of 
the CoRoT mission, through multicolor broad band photometric 
observations by CoRoT-Exoplanet team. The Exo-Dat catalog 
consists in 10.6 millions of stars over an area of 220 deg2, ob- 
served with the Wide Field Camera (WFC) at the 2.5m Isaac 
Newton Telescope (INT) at Roque Muchachos Observatory on 
La Palma, Canary Islands. Harris B and V filters and r' and i' 
ones from the Sloan-Gunn system, were used (the limiting mag- 
nitude is about r' - 20.0). The spectral type, luminosity class 
and the mean reddening were obtained by fitting models atmo- 
sphere to the observed fluxes (Deleuil et al. 2008 and M. Deleuil 
private communication). All stars classified as BO-AO (IV-V), 
with low contaminatiorQ, were systematically included in our 
programme, ensuring about 2 HgMn stars per set of targets. This 
programme was part of a larger one involving other teams con- 
cerned by this survejQ. 



^ CoRoT was launched on 2006 December 27. 

^ Only stars with low contamination are eligible for CoRoT obser- 
vations. Contamination (by surrounding sources) is determined in Exo- 
Dat through criteria established for planets detection. 

074.D-0193(A), Looking deep in tlie CoRoT eyes, G. Alecian (PI), 
A.-M. Hubert, M. Deleuil, C. Neiner, M. Floquet, C. Martayan. 



3. Ground data 

3. 1 . VLT-FLAMES/GIRAFFE observations 

Observations were done during the period 74A (0.5 night, 30-3 1 
January 2005), in the framework of the Paris Observatory GTO. 
In the present work, we use the data obtained in MEDUSA mode 
with setup LR02 (396.4-456.7nm, R=6400). Observational con- 
ditions were rather good with an average seeing of 0.6"-0.9". 
Two exposures of 30 minutes were done for each field and an 
average S/N of about 70 was reached for brightest targets. The 
spectra we use in this work were obtained after data reduction 
performed with the dedicated software GIRBLDRS developed 
at the Geneva ObservatorjB Several tasks of the IRAF package 
for extraction, calibration, and sky correction of the spectra were 
also used. 

3.2. Looking for HgMn stars 

About 240 stars were observed in this run, and about 60 of 
them were classified in Exo-Dat as BO-AO (IV-V) stars. To 
find HgMn candidates among these 60 stars, we searched for 
the characteristic pattern formed by two Mnll lines (420.53nm 
and 460.637nm). This pattern is clearly visible in the spectra of 
HD 175640 (a high S/N UVES spectrum of this star was kindly 
communicated by S. Hubrig, see Castelli & Hubrig 2004), and 
remains well recognizable after degrading the S/N to 100. After 
a close examination of the spectra of these 60 stars, one star 
(SI) was found to show the Mnll lines pattern very similarly 
to HD 175640 (degraded spectrum), and two other stars (S2 and 
S3) were found with the same pattern, but less clearly. The char- 
acteristics of these three stars are shown in Table[T] their number 
in USNO-A2 catalog, the visual magnitude from Exo-Dat, and 
the fundamental parameters discussed in Section 13.31 A rough 
estimate of various lines strength for usual abundances in HgMn 
stars was done using the COSSAM code (Stift, 2000), and we 
checked the consistency of our selection of HgMn candidates by 
considering the spectra at other wavelengths. These three stars 
were then proposed to be included in the target allocated to the 
APs (Weiss, 2006). 

The detailed abundance analysis of these three stars, in the 
next subsections, have confirmed their identification as good 
HgMn candidates. 

3.3. Fundamental parameters 

The effective temperatures (Tetf) and surface gravities (log g) of 
the stars were derived by adjusting a grid of theoretical Balmer 
lines to the observed one. With a lack in photometric data, fun- 
damental parameters can only be determined spectroscopically. 
By varying simultaneously Teft and log g, we select the pro- 
file that fits best Balmer lines. LTE model atmospheres were 
calculated using Kurucz's ATLAS 12 code (Kurucz, 2005), as- 
suming a plane parallel geometry, a gas in hydrostatic and ra- 
diative equilibrium and a depth-independent microturbulence. 
Abundances can be adjusted individually in ATLAS 12 which 
uses the Opacity Sampling technique for line opacity. The in- 
crement in Teff and log g were 500 K and 0.5 dex respectively. 
We have used SYNSPEC code (Hubeny & Lanz 1992) to com- 
pute the synthetic Balmer profiles in the wavelength range of 
GIRAFFE spectra (^3950-4550 A). H^ and H^ were generated 
using the line broadening functions of Vidal et al. (1973). The 



' http://girbldrs.sourceforge.net 
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Table 1. Programme stars and adopted effective temperatures, surface gravities and rotational velocities for our sample stars 



Star- 


USNO-A2 




T-eff (K) 


log g (dex) 


V'e sin ( (km.s ') 


SI 


0825-03036752 


12.05 


13500+500 


4.00+0.5 


35 + 5 


S2 


0825-03028353 


13.01 


12500+500 


4.50+0.5 


50 + 8 


S3 


0825-02993210 


13.18 


11750+500 


4.00+0.5 


50 + 8 




Fig. 1. The effect of changing Tefj by 500K and log g by 0.50 dex on the spectrum around H^. The adopted temperature and gravity 
for SI are 13500 K and 4.0 dex 



solar abundances used in SYNSPEC and ATLAS 12 are those 
from Grevesse & Sauval (1998). These Balmer lines are com- 
puted using a non solar ATLAS 12 model. After that we select 
a model atmosphere, we determine the abundances of the ele- 
ments and these abundances are included in the new computed 
model atmosphere for the same star and the final models have 
the abundances displayed in Table [3] 

For stars hotter than 8000 K, Balmer lines are more sensitive to 
gravity than to temperature. Figure [TJa shows the effects of in- 
creasing and decreasing the effective temperature by about 500 
K around 13500 K on the profile of H^. This variation of +500 
K corresponds to the estimated error on the effective tempera- 
ture. In the same manner, we have tested the effects of changing 
the surface gravity by about +0.50 dex which is the estimation 
of the error bar on log g . These tests were done for the three 
stars and for both Hy and H5 profiles. For stars earlier than mid- 
A type, hydrogen lines are relatively free from metal lines and 
thus an appropriate normalization to the continuum can be done. 
The adopted fundamental parameters are those of the model that 
best fit the overall shapes of Balmer lines. The adopted effec- 
tive temperatures, surface gravities and rotational velocities are 
collected in Table □ for SI, S2 and S3. 



3.4. Abundance analysis 
3.4.1. Input data 

Once the fundamental parameters were fixed, the same 
ATLAS 12 model atmospheres were used as an input in the syn- 
thetic spectra calculations. We have compiled a line list in the 
optical domain with the most accurate atomic parameters. For 
Mgll, Sill, Till, CrII, Fell, SrII and YII, we have selected lines 
(Table|2]i from the atomic transitions list of Gebran et al.'s (2008) 
online table after we have checked that these lines are good 
enough for abundances determination in case of late B-type stars 
(ie. no blending and not so weak). These authors have used 
Kurucz's gfall.daQ as an initial line list, and they have modi- 
fied the atomic parameters of the investigated chemical elements 
(wavelengths, excitation potential, oscillator strength and damp- 
ing constants) by checking for more recent laboratory determi- 
nations (as in NIST and VALD databases). For Hel, CII, PII, 
Mnll and Hgll, we have used Kurucz's data list. The helium 
abundances were derived using two lines of Hel at /14143.8 and 
/14471.5 using the broadening table of Barnar et al. (1969). We 
have also included data for hyperfine splitting for the selected 
transitions when relevant, in particular for Mnll. The spectral 
resolution of the spectra and the effect of rotational broadening 
prevent us from detecting the signature of hyperfine splitting. 



http://kurucz.harvard.edu/LINELISTS/GFALL/ 
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Mercury abundance was derived using the line of Hgll at 398.39 
nm with the oscillator strength from Sansonetti & Reader (2001). 

3.4.2. The method 

Our main interest in the abundance determination is to detect 
HgMn features because we cannot derive accurate chemical 
abundances using low resolution spectroscopy and low signal- 
to-noise spectra. We have derived individual chemical abun- 
dance of He, C, Mg, Si, P, Ti, Cr, Mn, Fe, Sr and Y by iter- 
atively adjusting synthetic spectra to the observed normalized 
spectra. The synthetic spectra were computed using Takeda's 
(1995) code. Takeda's procedure minimizes iteratively the chi- 
square statistic between the synthetic spectrum and the observed 
one (See Takeda's paper for a full description of the method). 
Apparent rotational velocity (ve sin ;) and microturbulent veloc- 
ity (^,) were derived using the strong Mgll triplet at 4480 A 
and a set of neighboring unblended weak and moderately strong 
Fell lines between 4490 and 4530 A as described in sec 3.2.1 
of Gebran et al. 2008. Allowing small variations around solar 
abundances of Mg and Fe, we iteratively fitted these lines leav- 
ing v'e sin /and ^,as free parameters. The weak Fe II lines are sen- 
sitive to rotational velocity and not microturbulent velocity, the 
moderately strong Fell lines are sensitive to the microturbulent 
velocity changes. The Mgll lines are sensitive to both Ve sin /and 
Each of the Mg and Fe lines yielded a set of values for log e, 
Ve sin /and which were in good agreement. 

Once Ve sin / and ^, are fixed, abundances have been derived 
for each line of each chemical element in our sample stars. With 
a resolving power of about 6500, errors on rotational velocities 
are more or less of the order of ~15%. Microturbulent velocities 
were found to be very small, which is common in HgMn stars 
(see Smith & Dworetsky 1993), thus we decided to fix them to 
zero for SI, S2 and S3. The found apparent rotational velocities 
(ve sin /) are displayed in column 7 of Table [T] 

The derived chemical abundances relative to the surQ for 
each lines used for this analysis are shown in Table |2l The fi- 
nal abundances for each star, and their respective errors are dis- 
played in Table [3] The abundance [X/H] does not include the 
rescaling of the H abundance in order to account for the He defi- 
ciency. We have assumed that the abundances derived from sev- 
eral lines for a chemical element, follow a Gaussian distribution. 
The errors on the element abundances are thus standard devi- 
ations. If the abundance is derived from only one line, we as- 
sign to it an error of 0.20 dex. We have also tested the effect of 
changing the atmospheric parameters on the derived abundance. 
Changing the effective temperature by an amount of + 500 K af- 
fect the mean derived abundances by an amount of ~0. 1-0.3 dex 
depending on the element and of course the HgMn trends are not 
affected. Figure |2] displays an example of the iterative fitting of 
the Mnll line at ~4519 A for SI. The synthetic spectrum with 
solar abundances is in dotted line and the manganese enhanced 
spectrum is in dashed line. In case of S 1 , manganese is about 300 
times higher than solar. We have displayed in Fig[3]the overall 
fit of the three star between 4000 and 4400 A. In this region, we 
have the two Balmer lines (Hy and Hg) and a large number of 
Mnll lines. In FiglH we show some examples of line fitting for 
SI, S2 and S3 and for several chemical elements as He, Hg, Fe 
and Mg. Due to the low signal-to-noise, some of the weak lines 
couldn't be analyzed as the titanium in SI or nickel in S2 and 



([§] = log(f )* - log(§)o) 



Table 3. Abundances (relative to the sun) table for SI, S2 and 
S3. 



otar 


C 1 


CO 

oZ 




Hel 


-1.53 


-1.51 


-1.17 


CTHe 


u.uy 


n ^ 1 


U.ZU 


CII 


-0.57 


-1.30 


-0.68 


(Tc 


U.zU 


n on 
U.ZU 


n on 
U.ZU 


Mgll 


-0.39 


-0.50 


-0.45 


O" Mg 


n 1 ^ 
U. i J 


n TO 


U.U / 


Sill 


0.35 


-0.61 


0.02 




U.zU 


n 1 n 


n 1 Q 
U.15 


PII 


2.16 


- 


- 


cr p 


U. iz 






Till 




1.23 


0.43 


C Ti 




n n^ 
U.Uj 


n 1 n 

u. iy 


CrII 


1.10 


-0.121 


0.48 


o"c> 


U.zU 


n on 
U.ZU 


n on 
U.ZU 


Mnll 


2.51 


1.83 


1.42 


O" Mn 


n i'2 
yj.Zj 


n 

U.Zj 


U.4 / 


Fell 


0.40 


-0.10 


-0.08 


cr ff 


n ]A 

U. 


U. 10 


n 1 ^ 
yj.ij 


Nill 


-0.75 








0.20 






SrII 


0.92 




2.44 




0.20 




0.20 


YII 




3.24 


3.65 






0.20 


0.20 


Hgll 


6.61 


6.11 


4.97 


O-Hs 


0.20 


0.20 


0.20 



1.05 



0.9 





1 1 1 1 1 1 1 

Observed spectrum 




II 




Solar SS 








Mn -enhanced SS 








\ 

A ^ ■■. / 


ij 


l\ / 




\\ /\'. / 
\\ / V. / 
A / V / 

\ \ 




\ I 
\ 1 

n / 
A / 
A / 




Mnll Fell 






\y 

Fell 
II 



451.8 4520 4522 4524 



A (A) 

Fig. 2. Adjustment of synthetic spectra on the observed one for 
SI. The dotted lines represent a synthetic spectra with solar 
abundance for all the elements. The dashed lines represent a syn- 
thetic spectra with an enhancement of manganese by about 300 
times the solar value. 



S3. It is clear that all these 3 late-type B stars present HgMn 
features. Large overabundances are found for P, Mn, Sr, Y and 
Hg and underabundances in He and Ni. It also appears that the 
hotter HgMn star have the stronger helium deficiency and man- 
ganese enhancement. Figure |5] displays the abundances pattern 
of SI, S2 and S3. One should be aware that the present analysis 
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Table 2. Detailed abundances (in log (Nx/Nh) + 12 scale) for each line of each chemical element. 



5 



Ion 


A (A) 




SI 


S2 


S3 


Hel 


4143.800 


-1.19 


9.486 


9.177 


9.500 


Hel 


4471.500 


0.05 


9.456 


9.800 


9.450 


CII 


4267.261 


-0.58 


7.906 


7.180 


7.804 


Mgll 


4390.510 


-1.70 


6.975 


7.359 


7.122 


Mgll 


4427.994 


-1.20 


7.270 


- 


7.002 


Mgll 


4481.10 


0.73, -0.57, 0.57 


7.201 


6.720 


7.154 


Sill 


4128.054 


0.30 


- 


7.090 


7.704 


Sill 


4130.894 


0.46 


7.857 


6.702 


7.350 


PII 


4420.712 


-0.48 


7.608 


- 


- 


PII 


4452.472 


-0.19 


7.420 


- 


- 


PII 


4466.140 


-0.28 


7.437 


- 


- 


PII 


4468.000 


-0.21 


7.760 


- 


- 


PII 


4475.270 


0.30 


7.670 


- 


- 


PII 


4483.693 


-0.43 


7.502 


- 


- 


Till 


4163.644 


-0.13 


- 


6.130 


5.880 


Till 


4287.873 


-1.79 


- 


6.206 


5.340 


Till 


4300.042 


-0.44 


- 


6.207 


5.246 


Till 


4386.844 


-0.96 


- 


- 


5.293 


Till 


4394.059 


-1.78 


- 


6.290 


5.450 


Till 


4395.051 


-0.54 


- 


- 


5.270 


Till 


4399.772 


-1.19 


- 


6.280 


- 


Till 


4417.714 


-1.19 


- 


6.180 


- 


Till 


4443.801 


-0.72 


- 


6.170 


- 


CrII 


4616.629 


-1.29 


6.740 


5.510 


6.113 


Mnll 


4085.390 


-2.51 


7.800 


- 


- 


Mnll 


4252.963 


-1.13 


7.820 


7.180 


- 


Mnll 


4259.200 


-1.58 


- 


7.089 


7.440 


Mnll 


4292.237 


-2.22 


7.310 


- 


- 


Mnll 


4275.884 


-1.91 


8.139 


- 


- 


Mnll 


4278.614 


-2.51 


- 


- 


7.080 


Mnll 


4308.158 


-1.72 


- 


7.68 


6.587 


Mnll 


4310.692 


-0.15 


7.677 


- 


- 


Mnll 


4325.041 


-2.29 


8.200 


- 


6.036 


Mnll 


4356.621 


-2.02 


7.809 


7.246 


- 


Mnll 


4363.255 


-1.90 


7.890 


6.980 


- 


Mnll 


4478.637 


-0.95 


8.040 


- 


6.708 


Mnll 


4503.201 


-2.16 


7.880 


- 


- 


Mnll 


4518.965 


-1.32 


7.850 


6.900 


- 


Fell 


4273.326 


-3.25 


7.866 


7.260 


7.496 


Fell 


4296.570 


-3.01 


- 


7.380 


- 


Fell 


4416.830 


-2.60 


7.600 


- 


- 


Fell 


4491.405 


-2.70 


7.887 


7.106 


7.123 


Fell 


4508.288 


-2.21 


7.790 


7.432 


7.399 


reli 


4515. 3jy 


-/.4o 


/.56U 


T C 1 A 

l.jlU 


n lACt 

/.34y 


Fell 


4520.224 


-2.60 




7.150 


7.346 


Fell 


4522.634 


-2.03 


7.760 


7.520 




Fell 


4541.524 


-3.05 


7.950 






Nill 


4067.031 


-1.83 


5.463 






SrII 


4215.520 


-0.17 


3.850 




5.367 


YII 


4374.935 


0.16 




5.440 


5.853 


Hgll 


3983.90 


-1.51 


7.700 


7.200 


6.065 



is based on low resolution and low signal-to-noise spectra, but 
the important result is that the HgMn pattern is clearly displayed 
in all three stars. 



4. CoRoT Observations 

Stars SI, S2 and S3 were observed by CoRoT during the first 
long run towards the galactic anticenter (LRaOl). This long 
run lasted from the 24th of October 2007 to the 3rd of March 
2008, i.e. about 130 days. All the results presented here are 



based on the N2 data provided by the CoRoT data centefl re- 
leased on October 29th, 2008. We have proceeded to the addi- 
tional treatments of the data (see Sec 14. It using IgorPro v6 (from 
Wavemetrics^ and Period04 vl.0.4 (Lenz & Breger, 2005) for 
the analysis of frequencies. 

In this section, for the sake of conciseness, we will present 
in details the treatment of data for the star S 1 only. 



http://idoc-corot.ias.u-psud.fr/index.jsp 
' http://www.wavemetrics.com/ 
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4.1. Data preparation 

A detailed presentation of the CoRoT data is given by Auvergne 
et al. 2009. For each of our three stars, CoRoT provided chro- 
matic lightcurves, in the oversampled mode (32 sec.|3 These 



Chromatic and oversampled data are available only for a limited 
number of targets of the exoplanet channel. 



chromatic lightcurves are not provided to give color informa- 
tions on targets, because the separation between the colors 
(which is realised by the ground based data processing) is ar- 
bitrary. The color separation was implemented to help in transit 
detection of exoplanets, but does not give workable information 
on star's properties, at least at the present stage of data treatment. 
The colored fluxes are strongly sensitive to depointing fluctua- 
tions and consequently the jitter corrections on colors are not 
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Fig. 5. Abundances pattern for our sample stars 

very efficient. The white Ught flux is less sensitive to jitter be- 
cause on the aperture edges the PSF slope is small. So it is rec- 
ommended to work first for detection on white lightcurves. 

The original set of N2 lightcurves (red, green, and blue) 
for SI consists in 338729 points of measurement (electrons vs. 
CoRoT- Julian days), corrected for main instrumental effects. We 
have first selected the points flagged as "valid" {status = 0) by 
the N2 pipe-line. Other points are generally affected by the so- 
called SAA (South Atlantic Anomaly) and/or other instrumental 
effects. However, some abnormal points survive this selection 
and need to be suppressed manually. At the end of this selec- 
tion, one has 300029 points of measurement (over 131 days). 
Most of them (98%) are spaced by 32 sec, but this whole se- 
lection process introduces some gaps, often smaller than about 
1100 sec due to SAA. The SAA is crossed by the satellite 8 
time per day, and lasts about 10 mn. There are also gaps due to 
various interruptions (normal and abnormal) of data acquisition. 
The largest gap is 13.7 hours from JD=2861.23 due to missing 
data. The next step consisted in correcting the monotonic long- 
term variation of each curve by a parabolic fitting. This variation 
is partly due to a pointing shift, and partly due to instrumental 
drift. After this correction we ended with three lightcurves sum- 
marized by the following values (average level, standard devi- 
ation): Red (391.373 kE, cr = 1.293 kE), Green (107.826 kE, 
o- = 0.435 kE), Blue (242.617 kE, cr = 1.056 kE). Finafly, 
for analysis of frequencies (Sect. |4.2| i, we have built the white 
lightcurve by simply summing the three chromatic lightcurves, 
and we achieved a binning of 512 sec to reduce the noise (here- 
after, all the presented plots and results are based on data with a 
binning of 512 sec). A similar treatment was done for stars S2 
and S3. 

4.2. Data analysis and results 
4.2.1. StarSI 

The white lightcurve (points) is shown in Fig|6] The solid line 
is the fit obtained with Period04, using the frequencies listed in 
Table |4] This list of frequencies corresponds to frequencies with 
the highest amplitudes, from which the well-known CoRoT- 
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Fig. 7. Full S 1 spectrum (amplitude vs. cycles per day) calculate 
with Period04 from the white lightcurve. 
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Fig. 8. Details of SI spectrum. Zoom in Fig|7] with a linear 
scale. 

Table 4. List of frequencies for star S 1 (from Period04), used to 
calculate the fit of Fig|6] 



N 


c/d 




Ampl. 


Phase 


Period (days) 


1 


0.0114 


0.132 


413.9 


0.7943 


88.007 


2 


0.0146 


0.169 


106.3 


0.8962 


68.412 


3 


0.0251 


0.290 


174.6 


0.0849 


39.880 


4 


0.0436 


0.505 


74.1 


0.7371 


22.913 


5 


0.0877 


1.015 


73.5 


0.4033 


11.401 


6(P) 


0.2324 


2.689 


552.9 


0.9211 


4.304 


7 (P/2) 


0.4648 


5.379 


212.7 


0.6286 


2.152 


8 (P/3) 


0.6982 


8.081 


37.1 


0.5796 


1.432 



orbital ones and its harmonics were removed. The full spec- 
trum is shown in Fig|2] The CoRoT orbital frequency at 13.97 
c/d (161.7//Hz) and its harmonics are clearly visible. The low 
frequencies around 0.01 c/d are related to the long-term varia- 
tion noticeable in Fig|6l We cannot interpret safely these long- 
term variations which have periods comparable to the run du- 
ration. They do not correspond to known instrumental effect. 
Fig. [8] presents the detailed spectrum around 0.23 c/d. The fea- 
ture around the peaks are due to the spectral window. In our opin- 
ion, the period P ^ 4.3 days, with two harmonics (P/2, P/3) is 
an intrinsic variation of star SI and its amplitude (Am/m) can be 
estimated to be about 1 .6 mmag. There is not any known CoRoT- 
periods around P, P/2 and P/3. It is easily visible in Fig|6](the 
"high" frequency variation). 

4.2.2. Star S2 

The lightcurves of star S2 are much noisier than those of star S 1 . 
In addition, the blue lightcurve presents some strong disconti- 
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Fig. 9. Full S2 spectrum (amplitude vs. cycles per day) calcu- Fig. 10. Full S3 spectrum (amplitude vs. cycles per day) calcu- 
lated with Period04 from the white lightcurve. lated with Period04 from the white lightcurve. 



nuities and perturbations from the date JD x 2931.36. We have 
then selected only points between JD - 2853.036 and 2931.35. 
These data lead to the rather noisy spectrum shown in Fig. |9] 
Again, the CoRoT orbital frequency (and harmonics) dominates 
at 13.97 c/d (and above). The only other significant peak is for 
0.337 + 0.01 c/d, very close to a period of 3 days, which may 
be related to the 1 day period often present in CoRoT data. We 
cannot find confidently from these rather poor data frequencies 
intrinsic to the star S2. 



4.2.3. Star S3 

As in the previous case, lightcurves of star S3 present several dis- 
continuities and perturbations, but not exactly at the same dates. 
However, compared to star S2, the lightcurves are less noisy 
and present interesting features. We have considered the data 
between JD = 2853.036 and 2925.933 (73 days) and we show 
the corresponding full spectrum in Fig. [TOl the details at low 
frequency range (< 3 c/d, linear scale) are in Fig.[TT] There are 
clearly two peaks above the noise, the first one at 0.395 ±0.01 c/d 
(P * 2.53 days) and the second one at 0.789 + 0.01 c/d which is 
possibly the first harmonic (P/2). There is a peak around 0. 19 c/d. 



P = 2.53 days 



0.45 day 




0.0 0.2 0.4 0.6 



1.8 2.0 2.2 2.4 



Fig. 11. 

scale. 



Details of S3 spectrum. Zoom in FiglTOl with a linear 



but as one can see in better detail in Fig.[TOl it does not emerge 
from the surrounding patterns as clearly as the previous peaks. 
Another isolated peak emerges around 2.24 c/d (0.446 day), but 
it seems less significant than the ones at 0.395 and 0.789. 
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5. Theoretical models for HgMn stars 

The models for this paper are calculated with atomic diffusion 
that is known from first principles. In radiative zones this implies 
taking into account atomic diffusion including gravitational set- 
tling, thermal diffusion and radiative accelerations, in addition to 
the purely diffusive term. The interactions between different dif- 
fusing species are also taken into account. Convection and semi- 
convection are modeled as diffusion processes as described in 
Richer et al. (2000) and Richard et al. (2001). The detailed treat- 
ment of atomic diffusion is described in Turcotte et al. (1998). 
In these models the Rosseland opacity and radiative accelera- 
tions are computed at each time step in each layer for the exact 
local chemical composition using monochromatic opacities of 
the OPAL group Iglesias & Rogers (1996). The radiative accel- 
eration for each element is from Richer et al. (1998) with cor- 
rection for redistribution of momentum to electrons and others 
ions from Gonzalez et al. (1995) and LeBlanc et al. (2000). All 
models have an initially homogeneous solar composition. 

For this preliminary study we include mixing in all the mod- 
els in an ad-hoc parametric way as in Richer et al. (2000) where 
the parameters specifying effective transport coefficients are de- 
fined. All our models were computed with a effective transport 
coefficient, Dj, 50 times larger than the He atomic diffusion co- 
efficient at log r = 5.3 and varying as p^^ deeper in the star (for 
temperature below logT - 5.3 we assume a complete mixing). 
The choice of the layer with log T = 5.3 as lower boundary of 
the mixed zone is justified by the fact that the iron accumula- 
tion (due to diffusion process) around that layers, leads to the 
appearance of a convection zone. The existence of HgMn phe- 
nomenon shows that full mixing cannot extend deeper, otherwise 
superficial abundance anomalies would not appear According 
to the diffusion model, this lack of deeper mixing may be due 
to the slow rotation of HgMn stars. To better model HgMn 
stars one would like to compute atomic diffusion in higher lay- 
ers (with lower temperature, including the atmosphere). But the 
monochromatic opacities available to us are not sufficiently ac- 
curate to compute reliable radiative accelerations at low temper- 
ature (LeBlanc et al., 2000; Richer et al., 1998) and also, mov- 
ing the boundary of the fully mixed zone to lower temperature 
leads to numerical difficulties in the present version of our code. 
Therefore, in our models, mixing impedes atomic diffusion in 
layers with log T < 5.3. Notice that we have not adjusted the 
mixing depth to fit the surface abundances of our observed stars, 
because these anomalies are expected to be due to atomic dif- 
fusion in the atmosphere. For layers deeper in the star, where 
logT > 5.3, the effective transport coefficient is adjusted to 
rapidly decrease the efficiency of the mixing so in these layers 
the effect of atomic diffusion leads to chemical stratification and 
to an accumulation of iron which affect the opacities and the pul- 
sational properties of the models. 

We have computed models of 3.OM0 and 4.QMq which are 
in the mass range of SI and S3 stars (see Figure [T2]|. As shown 
in Figure[T3]these models have a important bump in the opacity 
profile around log T - 5.3, at this temperature radiative accelera- 
tions lead to an iron overabundance. In this region, see Figure[T4l 
oscillation could be driven by /f-mechanism as -^{kt + f^) > 0, 
where /cj- = {dlnK/dlnT)p and^p = (dlnK/dlnpjr- 

The seismic properties of the models were studied with the 
help of the nadrot code (Dziembowski, 1977). Both the 3.OM0 
and the A.OMq models show positive growth rates for non radial 
gravity modes, and all p-modes are stable. In the A-.QMq model, 
the periods can be as short as 16 hours and as long as 2 days for 
young models, whereas no mode under 33 hours is overstable 




Fig. 12. Evolution of models of 3.OM0 and 4.OM0 in the logg - 
Teif plan. SI, S2 and S3 are shown with eiTor bars. 
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log T (K) 

Fig. 13. Rosseland mean opacity versus temperature in the 
4.OM0 model at the age of 100 Myr. 




5 5.5 6 

log T (K) 

Fig. 14. Logarithmic opacity derivatives vs temperature for the 
4.OM0 model at the same age as in Figure [13] 



at the end of the main-sequence. The upper limit for the period, 
when the 4.OM0 leave the main-sequence, is hard to define as 
the number of modes identified by the pulsation code is limited, 
but we find modes with period as long as 4.4 days to be over- 
stable while modes with period above 6.5 days are stable. The 
3.OM0 model have a similar range in periods for young models 
but show a smaller shift to longer periods during evolution with 
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no overstable modes with period longer than 2.5 days until all 
modes become stable around 250 Myr. 



6. Discussion 

The variations we find for stars SI and S3 can be reasonably 
considered as intrinsic to these stars. 

The VLT/GIRAFFE spectra used in this work are medium 
resolution spectra, then the v.sini (Table [TJ we obtained are not 
very accurate and can be rather considered as upper limits of 
the values one could expect to determine from higher resolution 
spectra. Since v.sini fixes a lower limit of equatorial velocity, 
it gives an indication about the upper limit of rotational period. 
The upper limits of rotational period which can be inferred from 
Table [1] appear to be in the same ranges as the periods we find 
for SI and S3. Therefore, one cannot exclude at this stage that 
the variations found in our lightcurves are due to rotational mod- 
ulations (see for instance, Lanza et al. 2008). 

6.1. Rotational modulation hypothesis 

According to the discussion by Clarke (2003), surface inhomo- 
geneities like spots or patches can produce photometric varia- 
tions, but only the fundamental and the second harmonic could 
be present and not the third harmonic as we observe in SI. 
However, it is difficult to assert that high degree harmonics can- 
not be found in Fourier analysis in case of complex superfi- 
cial structures. Magnetic Ap stars are known to have large scale 
and organized strong magnetic fields, with abundance patches of 
many metals. These anomalies are qualitatively well understood 
in the framework of atomic diffusion but the detailed quantitative 
modelling of that process is difficult, and sometimes strongly 
suggests that diffusion has to work with other processes like in- 
homogeneous winds (Babel 1992). One notes also that a com- 
plete discussion should also consider light-induced drift and am- 
bipolar diffusion. Such superficial structures have been studied 
from several decades from line shape variations (for recent stud- 
ies, see for instance observations by Kochukhov et al. 2004, and 
modelling by Alecian & Stift 2007), and one would not be sur- 
prised if such complex structures cause photometric variations 
(in addition to those due to the roAp phenomenon for cool mag- 
netic Ap stars). If our stars were not HgMn, but magnetic Ap 
stars, one would consider attentively rotational modulations to 
explain lightcurve variations. 

Considering that the 3 stars discussed in this work are likely 
to belong to the HgMn group, is the hypothesis of abundance 
patches or spots relevant in explaining the photometric varia- 
tions we observe? Observational facts are that there is not yet 
any confirmed detection of magnetic fields (Shorlin et al. 2002, 
Wade et al. 2006) nor any detection of photometric variations 
(Adelman et al. 2002) in HgMn stars from ground based obser- 
vations. There are some possible cases of few hundred Gauss 
detections (Hubrig et al 2006a), but they need to be confirmed, 
and they are not typical for the HgMn group. Following these 
constraints, standard diffusion modelling of HgMn stars consists 
in estimating the diffusion velocities and the resulting vertical 
stratifications of elements (equilibrium approximation), assum- 
ing stable atmospheres without magnetic fields. This modelling 
produced encouraging results by explaining main trends of metal 
overabundances in these stars (Smith & Dworetsky, 1993), but 
cannot yet reproduce the large variety of detailed observations 
of individual stars. An important reason of the difficulty in mod- 



elling individual stars is that diffusion is a slovvFI (and then frag- 
ile) time-dependent process which, at a given time, has built up 
abundance stratifications which are different from those com- 
puted assuming equilibrium (the finale stage of an ideal unper- 
turbed time-dependent process). 

Because, in the standard model, atomic diffusion is supposed 
to work alone, any perturbation, as it could exist in a real HgMn 
star, can modify the picture of the stratification build-up pro- 
cess. For instance, effects induced by binarity, slow matter cir- 
culation, weak turbulence or wind can exist in HgMn stars. If 
so, they could be undetectable but strong enough to affect the 
stratification process. Among possible theoretical scenarios, one 
can consider the existence of weak, unorganized magnetic fields 
(hardly detectable). In magnetic Ap stars, where magnetic fields 
strengths can be larger than several kG, the component of ion's 
diffusion velocities orthogonal to the field lines is significantly 
reduced above logr » and vanishes as matter density de- 
creases. This leads to the formation of non-uniform (horizon- 
tally) clouds of metals, anchored in the usual line-forming re- 
gion (just above logr w 0). In HgMn stars, if weak fields exist, 
the effect on diffusion velocities will be the same as in magnetic 
stars, but much higher in the atmosphere (for instance, above 
log T ^ -4.0 for fields < 1 kG, see Alecian & Stift 2006). This is 
due to the fact that the diffusion velocity of an ion in a magnetic 
field is slowed down by a factor which varies more or less as the 
square of the proton number density (see the detailed formulae 
in Sec. 4.1 of Alecian & Stift 2006). So, a weak magnetic field 
which will have no noticeable effect in usual line-forming re- 
gion could, however, force elements to remain bound to the star. 
In that case, one could expect that there may be high altitude ad- 
ditional clouds for some elements, forming abundance patches 
(related to the magnetic field orientation, possibly unorganized) 
and superimposed to the uniform stratifications predicted by the 
standard diffusion models. 

Recent works (Adelman et al. 2002, Kochukhov et al. 2005, 
Hubrig et al., 2006b) have detected in some HgMn stars varia- 
tions in some sparse line profiles and interpreted them in terms 
of abundance spots. As discussed just before, these observations 
are not necessarily in contradiction with the diffusion model. 
Moreover, one can note that these variations of line profiles are 
essentially found for elements having low cosmic abundances 
(Hg, Y, Pt, Zr), which bring an additional support to the diffu- 
sion model. Indeed, according to what one knows about radia- 
tive accelerations, these kind of elements are generally strongly 
pushed up by the radiation fielcQ and can more easily accumu- 
late than abundant elements in the upper atmosphere (in addition 
to deeper stratifications), where a weak magnetic field becomes 
able to impede atomic diffusion to expel that element toward the 
circumstellar medium. 

According to the published profile variations mentioned 
above (papers do not give any estimation of the photometric vari- 
ations that can be induced by these profile variations), one can 
hardly consider that, even if our stars would have same kind of 
line profile variations, they could produce the 1 .6 mmag photo- 
metric variations we observed in CoRoT lightcurves. Remember 
that the CoRoT lightcurves are integrated photon flux from about 
400nm to lOOOnm. To have regular photometric variations (as 
the ones we observed) from abundance spots, equivalent width 
of lines should have to vary more or less in phase. This is not ob- 



" Diffusion time scales in atmospheres extend from some years to 
centuries according to the considered element. 

Because photoabsorption by bound-bound transitions remain unsat- 
urated until huge overabundances are obtained. 
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served and this is not consistent with abundance patches which 
will be distributed differently according to elements. It does not 
mean that line profile variations could not coexist with photo- 
metric variations. We just say that these two types of variations 
are related to two different mechanisms. 

Our conclusion is that the variations observed recently in 
some sparse line profiles for some HgMn stars remain consistent 
with the diffusion model, but could hardly explain the photomet- 
ric variations we observed in CoRoT lightcurves. Nevertheless, 
even if our feeling is that these variations are not rotational mod- 
ulations, more precise measurements of v.sini are clearly needed 
to discriminate pulsations and rotational modulations, and check 
whether these stars are multiple or not. This is planned for the 
near future. 

6.2. Pulsation hypothesis 

In the hypothesis of pulsations, we have a significant 4.3 days pe- 
riod in SI, and 2.53 days period for S3. They are compatible with 
periods of non-radial g-modes predicted by models presented in 
Sec|5] According to these models and the study of Turcotte & 
Richard (2003), these pulsations should be of the same type as 
those of SPB stars (Dziembowski et al. 1993). However, one can- 
not consider these stars as belonging to the SPB class, since stars 
of this latter class have more or less solar abundances. In average 
SPB stars are hotter than HgMn stars, and have higher gravity. 
But, the boundary at the cool edge and at the low gravity edge 
of SPB's (De Cat, 2007) overlaps respectively the hot/large ones 
of HgMn group. In the same way, the periods we find for S 1 and 
S3 are little bit longer and with smaller amplitudes (less than 
or around 1.6 mmag) than what is generally found for SPB's, 
but here again these variations are compatible with those ob- 
served for some modes of SPB's. If pulsations were confirmed 
in the future in HgMn stars, differences could come from the 
lack of superficial convection zone which leads to strong metal 
stratifications in outer layers in HgMn stars(including the atmo- 
spheres). One can speculate that could cause stronger damping 
of those modes compared to SPB stars. 



7. Conclusions 

In this work, we first identified from multiobject spectroscopy 
with moderate resolution (using VLT-GIRAFFE spectrograph), 
3 faint stars presenting the characteristic features of HgMn types, 
and eligible for observation with the exo-CCD's of the CoRoT 
satellite. The analysis of their lightcurves, obtained during the 
run LRaOl of CoRoT, shows periodic variations for at least two 
of theses stars. The periods of these variations are compatible 
with theoretical predictions. This could be the first detection of 
pulsations in HgMn stars. However, further ground based high 
resolution spectroscopic observations are needed to confirm the 
HgMn nature of these stars and that these variations are pulsa- 
tions and not rotational modulations. 
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